Transverse multi-mode effects on the performance of photon-photon gates 
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The muhi-mode character of quantum fields imposes constraints on the implementation of high- 
fidelity quantum gates between individual photons. So far this has only been studied for the lon- 
gitudinal degree of freedom. Here we show that effects due to the transverse degrees of freedom 
significantly affect quantum gate performance. We also discuss potential solutions, in particular 
separating the two photons in the transverse direction. 
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I. INTRODUCTION 

Photons are attractive as carriers of quantum infor- 
mation because they propagate fast over long distances 
and interact weakly with their environment. Their utility 
for quantum information processing applications such as 
quantum repeaters [H or quantum computing [2] would 
be further enhanced if it was possible to efficiently im- 
plement two-qubit gates between individual photons. 
Such two-qubit gates can be implemented probabilisti- 
cally using just linear optics and photon detection [3|, 
but strong photon-photon interactions would allow much 
more direct and deterministic implementations. Several 
approaches to the implementation of interaction-based 
photon-photon gates have been proposed, including pro- 
posals based on Kerr non-linearities in fibers or crystals 
[J], on electromagnetically induced transparency (EIT) 
in atomic ensembles and on the interaction of both 
photons with an individual quantum system Q. See Refs. 
|7h11| for recent experimental progress. 

In the simplest case, an ideal controUed-phase gate 
performs the transformations |0)|0) |0) |0), |0) |1) — 
|0)|1),|1)|0) ^ |1)|0),|1)|1) ^ e^1l)|l), where |0) and 
|1) are zero- and one-photon states respectively. In the 
present context, the phase (p acquired by the state |1)|1) 
is due to the interaction of the two input photons. For 
quantum information processing applications it is desir- 
able to achieve (j) = n. In the first theoretical papers 
the main focus was on how to achieve interactions that 
are sufficiently strong to allow large phase shifts. The 
photonic pulses were typically idealized as single-mode. 
Later it was realized that the (longitudinal) multi-mode 
character of the pulses imposes important constraints on 
the implementation of high-fidelity quantum gates [l2| - 
[lOj . The phase shifts due to the interaction depend on 
the relative position of the two photons, and take differ- 
ent values over the pulses because the photons have to 
be described as extended wave packets rather than point 
particles Hq]. As a consequence, an initial product state 
of the two photons is mapped by the interaction onto 
an output state that exhibits unwanted entanglement in 
the photons' external degrees of freedom. For large phase 
shifts this leads to low fidelities for the simplest quantum 
gate proposals [isj . 

It has been suggested that this difficulty can be over- 



come by more sophisticated quantum gate designs where 
the two photons pass through each other. This can 
be achieved by trapping one [1, [11] or both [m] of the 
two photons, by having a counter-propagation geome- 
try [isl ITg} , or by considering two photons with different 
group velocities [Toj . 

Here we show that having the photons pass through 
each other is not sufficient on its own, due to the pres- 
ence of the transverse degrees of freedom. In short, the 
interaction-induced phases also depend on the relative 
transverse position of the two photons, which leads to 
fidelity limitations that cannot be mitigated by counter- 
propagation. In the following we describe these limita- 
tions in detail. We also discuss potential solutions, in 
particular separating the two wave packets in the trans- 
verse direction, which is possible for non-linearities based 
on long-range interactions. 



II. INTERACTING PULSE EVOLUTION 

The general picture for photon-photon gate is the in- 
teraction between two fields (x, t) and ^2 (x, t) in three 
spatial dimensions. The field operators at t = are 
defined as ^'i(x) ~ '^J2k^iM^^^^^ where V is the 
quantization volume. They might describe photons in 
a Kerr medium or polaritons in an EIT medium. The 
field operators satisfy the equal-time commutation rela- 
tion [«',(x,t),^](x',t)] = (5y(5(x-x'). This means that 
photon absorption is negligible (e.g. in the EIT case the 
spectra of the pulses are well inside the transparency win- 
dow), so the evolution of the interacting fields can be 
regarded as unitary. Making the standard slowly vary- 
ing envelope and paraxial approximations, we have the 
effective Hamiltonian ^] (ft = 1 is adopted hereafter) 



K 



E 



d^xi'\(yi,t){v-V, 
■' i 



^}*.(x,t) (1) 



to describe their free evolution. Here v is the group ve- 
locity in the positive or negative z direction, and k = ^ 
is the carrier wave vector, is the transverse Laplace 
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operator. The interaction for the two fields is 

V = J ^'^1 / rf'a;2{«'l(xi,t)^](x2,i)A(xi-X2) 

*,(x2,t)#,(xi,i)}, (2) 

where the terms of i = j and i ^ j 'm the sum correspond 
to self-phase and cross-phase modulation effect, respec- 
tively. Given the total Hamiltonian H = K + V, the 
equation of motion, idt^i{:x.,t) — [^i{x,t), H], for the 
counter-propagating fields read 

(757 - *«^)*2(X2, = -ld(x2, 0*2(X2, t) 

(3) 

with the interaction potential 

d(x„0= /dVA(x, -x')(/i(x',t)+/2(x',t)), (4) 



where /„(x) = *t (x)*„(x). 

A photon-photon gate is implemented by evolv- 
ing an input bi-photon state |$) — |l)i|l)2 = 
/ d3xi/i(xi)*l(xi) / d3a;2/2(x2)*t(x2)|0), where 
/i(x) — (0|^'i(x)|l) are the pulse profiles, under the 
unitary time evolution 



U{t) = Te-'Jo 



z /„' df'Hif) 



-zf^dt'K-tJ*dt'V-zC 



(5) 



where T denotes the time-ordering operation, and the op- 
erator is factorized using the Baker-Campbell-Hausdorff 
formula. The first factor in Eq. ([5]) describes the free evo- 
lution including pulse propagation and pulse diffraction. 
The second factor is from the interaction between pulses. 
The third factor contains all commutators between the 
exponents of the first two terms; they reflect the inter- 
play between pulse motion and pulse interaction, which 
generally changes the pulse profiles. 

The ideal output state under a gate operation would be 
g'^g-' /o d,t ^1$^^ with (j) being a homogeneous controlled 
phase, where we are taking into account the free evolution 
of the photons. The actual output state, however, will 
heUim) = / d^xi f d3a:2V'(xi,X2,t)^l(xi)*|(x2)|0) 
[23| , giving the two-particle wave function 
V'(xi,X2,t) = (0|^'i(xi,t)^'2(x2,i)|$), which is generally 
non-factorisable with respect to Xi and X2. The fidelity 
F and the controlled phase of a gate operation are 
determined via the overlap between the actual output 
state and the freely evolved state 1$/?) = e~^-f^o ''■*'^\^): 

d^Xid^a;2i/'o (xi, X2, t)V'(xi, X2, t), (6) 



where ?/'o(xi, X2, is the corresponding two-particle 
function for the freely evolved state | ^ji) . 



The field equations ([3]) allow one to obtain the evo- 
lution of V'(xi,X2,t) by multiplying 'I'2(x2,t) to the 
right of the first equation of (3) and >I'i(xi,i) to the 
left of the second; then, the product with (0| and |$) 
is taken on the addition of the equations, yielding the 
following linear equation for the two-particle function 
(0|*i(xi,t)*2(x2,t)|$): 

ot ozi 
+ (^-^;^)(0|#i(xi,t)*2(x2,t)|$) 



{iv 



V2 



V2 

■ll^^)(0|*i(Xi,t)*2(x2,t)|$) 



2fco 2fco 
= -^A(X2-Xl)(0|*l(xl,^)#2(x2,^)|$). (7) 

Here we have used ^'z(x)|0) = 0, *,(x, t)5'|"(x', t)|0) = 
V Eke*'''^''"'''^|0) = ^(^Hx - x')|0). With the center of 
mass coordinate, X = ^^2^^ ^^'^ x = Xi— X2, the derived 
equation for i/)(xi,X2,t) — i?(X, t)^(x, i) can be reduced 
to 

with i?(X, t) being trivially evolved under the diffraction 
term e"*^^T,x, 

The above linear equation greatly simplifies the de- 
termination of the evolution of interacting pulses. This 
simplification is possible because we are considering the 
interaction between two single photons (as opposed to 
multi-photon pulses). In the co-moving coordinate that 
eliminates the term 2vdz£,{x,t) in ([5]), the three factors 
in Eq. ([S]) are translated into e~^^. 



-iC' 



respectively, to evolve ^(x, 0). Now the third factor 
contains the exponentials of the commutators between 

fl^Ta; ^^'^ fi^jt) as follows: 



6[2fc^^-[^'2fc^^J> 
(9) 



These commutators in the exponentials are of order l/r, 
where I = vt is the medium length and r = ka^ is the 
Rayleigh length, with a the transverse size of the pulses 
at f = 0. It is not difficult to achieve Z/r <C 1 in practice, 
making the effects from the third factor insignificant. For 
simplicity we will not include the first and third factor 
in the two-particle functions derived below, but the first 
order contribution from the third factor will be consid- 
ered in our numerical calculations. The second factor 
g-i¥=(x,t)^ where <p(x,t) = J^dt'A{xT,z - 2v{t - t')), is 
the main concern of the present paper. We will explain 
its effects with two examples. 
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FIG. 1: This qualitative plot shows the two-particle wave 
function as a function of the transverse coordinates x\ and 

X2 of the two photons. The photons are distributed over a 
multitude of transverse positions. For a contact interaction 
the interaction induces a non-zero phase only if the two pho- 
tons happen to be at exactly the same position (the diagonal 
line xi — X2 in the figure), see Eq. (10). The probability of 
this occurring is infinitesimal. As a consequence, the effective 
output phase defined by Eq. (6) is always zero, independently 
of the strength of the interaction. 
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FIG. 2: Schematic setup for a photon-photon gate working 
with counter-propagating single photon pulses m a medium 

under EIT conditions. The pulses interact with each other 
through the dipole-dipole force between the Rydberg states 
\di). The inset shows the relevant energy levels of the atoms. 
The pulses collide head-on, but they have a transverse ex- 
tent a. This leads to a dependence of the interaction-induced 
phase for the two-particle wave function on the relative trans- 
verse position, resulting in a trade-off between the effective 
phase (j) of the two-photon operation and its fidelity F. Only 
very small phases are compatible with high fidelities. 

IV. DIPOLE-DIPOLE INTERACTION 



III. CONTACT POTENTIAL 



Our first example is a highly local interaction described 
by a delta function potential, A(xi — X2) — Vb^'^-'(xi — 
X2). This is a good model for non-linearities that are due 
to the interaction of both photons with the same atom in 
an atomic ensemble or crystal, or to short-range atomic 
collisions 0, S [Tsl - fisj . This interaction gives an output 
two-particle function 



V'(xi,X2,t) = /i(zi - vt,yLTs)h{z2 + vt,yLTa) 
X exp{i^[i7(z - 2vt) - H{z)]5^^\i^t)}. (10) 



where H{z) is the Heaviside step function. From Eq. 
(ITO)) one sees that the interaction-induced phase is non- 
zero only if the transverse coordinates of the two photons 
coincide, i.e. on a subset of configuration space x^ = 
that has measure zero, see Fig. 1. As a consequence, 
one has F = 1 and ^ = 0. In the case of an ideal three- 
dimensional delta function potential the effective output 
phase is exactly zero, no matter how strong the interac- 
tion between the two photons. This is closely related to 
the results of Refs. [ij, [l7| for the one-dimensional, but 
co-propagating case. One finds essentially equivalent re- 
sults for any interaction whose range is much shorter than 
the transverse size of the wave packets. Note that this re- 
sult is consistent with the non-zero conditional phase for 
photon-photon interactions obtained in Refs. [l^, [25j , 
where the evolution is non-unitary, as manifested by a 
different field operator commutator. 



Related difficulties arise also for more long-range in- 
teractions. This can be seen from our second example, 
which is motivated by Refs. [HI, It concerns the in- 
teraction between polaritons whose atomic component is 
in a highly excited Rydberg state in an external electric 
field, cf. Fig. [21 This induces a dipole-dipole interaction 
between the polaritons. 



A(xi - X2) = C(l - 3cos' ?9)/|xi - X2|^ 



(11) 



where C depends on the specific Rydberg states used 
and -i? is the angle between Xi — X2 and the external field 
(along which the electric dipoles of the Rydberg states are 
aligned). This is an attractive system because Rydberg 
states have large dipole moments, leading to potentially 
very strong interactions between the polaritons [26l - [28j . 

We consider the situation where the external field is 
perpendicular to the direction of motion. We assume the 
initial pulse profiles to be /i(xi) = ■!/'o(a^i)V'o(yi)V'o(-Zi) 
and /2(x2) = ■0o(a;2)V'o(y2)V'o(z2 - I), where ipnix) = 
[^^^]5i/„(f )e-^(f with i/„(f ) being the Hermite 
polynomials. The evolution according to Eq. gives 
the output two-particle wave function 

/i(^i - Z,XT4)/2(^2 + ;,XT,2)e-^^("^'"^''/"^ (12) 

The interaction-induced phase in the above is given by 



ip{z,yiT,l/v) 



C 1 , ^ 2zx| 



-{- 



2V X^ ^(z2+x2,)i 

(z - 2/)3 + 2(z - 20x|, 
((z-202+x|)i 



}. (13) 



By Eq. ^ , the conditional phase and the fidelity F 
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in this case are determined as follows: 



J x'j, 

^^^ {z' + lf + 2{z' + _ {z' -If + 2{z' 



((z' + O^+x'^) 



where x'^ = xi — Icz, X2 



((z'-/)2+x'^)S 



(14) 



X2 , and x' = x'j^ ^ Xg . In the 



calculation we have chosen I = A-Ka and a = 10 A. The 
results are shown in Fig. [2j Increasing the parameter 
C/{2va'^), which indicates the interaction strength, in- 
creases the effective output phase 0, but diminishes the 
output fidelity F. As a consequence, significant phase 
shifts are completely out of reach if one wants to achieve 
high fidelities. In the numerical calculations, we have in- 
cluded the first order correction due to the third factor 
of Eq. (I5|). Such correction comes from the commutators 

and [[-jp, lys], see Eq. (|9]). The commutators 
involving higher powers of if are shown to be vanishing. 
The exponentials of these commutators effect a modifi- 
cation of the intensity profile (x^) and a modification 
of the phase profile e~^'^^^'^^\ respectively. 

The main cause for the trade-off between cf) and F is the 
dependence of the interaction-induced phase (p{z, x^) on 
the transverse relative position x^. In fact, most of the 
behavior shown in Fig. [5] can be understood by setting 



the phase in the integrand of Eq. (ITi|) as 



C 1 



is quite accurate for Z 3> cr. 
mode mixing in the form 



which 



It gives rise to transverse 



V'o (a^i ) V'o (yi )i/'o (a;2 ) V'o (2/2 ) 



{^C„in1pmixi)lpn{yi)}{^dik1pi{x2)'ipkiy2)}, (15) 
m,n l^k 

leading to the deviation from the ideal output two- 
particle function e*'^Vo(2:i)V'o(yi)'0o(a;2)'0o(2/2)- 

The above analysis shows that the transverse mode 
mixing (or transverse mode entanglement) will develop 
even if the pulses are initially in a single transverse mode. 
This effect, which significantly affects the performance of 
photon-photon gates, is here described for the first time, 
to the best of our knowledge. Our analysis also clarifies 
that the pulse diffraction in the transverse direction has 
no direct impact on the performance of a photon-photon 
gate. It infiuences gate performance through its interplay 
with the interaction between pulses, i.e., by the third fac- 
tor in Eq. ([S])- Compared with the effect of the trans- 
verse mode mixing shown in Eq. (|15p . such diffraction- 
interaction interplay is insignificant in the regime con- 
sidered here, where the medium length I is much smaller 
than the Rayleigh length r. For example, for l/r — 0.2 
(as in our calculations) it induces corrections only at the 
few-percent level. 
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FIG. 3: Introducing a transverse separation between the two 
pulses greatly relaxes the trade-off between F and (f), such that 
a phase of order tt becomes compatible with high F. The di- 
mensionless separation R is defined as The price to pay is 
that the interaction strength has to be increased significantly 
in order to compensate for the transverse separation. 



V. POTENTIAL SOLUTIONS 

We have seen that the transverse multi-mode charac- 
ter of the quantum fields, which leads to a transverse 
relative position dependence of the interaction-induced 
phase shifts, has very significant consequences for the fi- 
delity and phase achievable in photon-photon gates. We 
will now discuss two potential solutions for this prob- 
lem. The first is applicable only to the case of long- 
range interactions. It consists in separating the paths 
of the two photons by a transverse distance D that is 
much greater than the transverse size a of the pulses, 
i.e., the initial profiles of the pulses will be, for exam- 
ple, /i(xi) = ipo{xi)i'o{yiHo{zi) and /2(x2) = tpo{x2 - 
D)ipoiy2)''pQ{z2 — I) ■ With increasing D one will approach 
a situation where the transverse degrees of freedom of the 
photons can effectively be treated as point-like. As a con- 
sequence, the effect studied here will diminish. This is 
shown in Fig. [3l We have again chosen a medium length 
I = Attu. Interaction-diffraction interplay effects are at or 
below the 10~^ level in this case because they depend on 
the gradients in the interaction-induced phase across the 
wave packets, which decrease with increasing transverse 
separation. 

The price to pay for the improved fidelities shown in 
Fig. [3] is that even for dipole-dipole interactions the 
interaction-induced phase decreases quickly as the trans- 
verse separation is increased. For example, achieving a 
conditional phase shift = tt with a fidelity F = 0.9 
requires R = — = 26. Comparing to Ref. [l^, this 
means that one could work with a principal number of 
the Rydberg state n ~ 75, a transverse wave packet size 
a = 7^m, and a group velocity v = 4 m/s. Achieving 
(j> = TT with a fidelity F = 0.99 requires R — 79, which is 
possible provided that, for example, n can be increased 
to 100, V reduced to 1 m/s, and a reduced to 5 /im. 
These requirements are realistic with current technology. 
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in particular Rydberg states with n = 79 were already 
used in the experiment of Ref. [1^ . 

The second potential solution, which is applicable both 
to short-range and long-range interactions, consists in im- 
posing strong transverse confinement. If the confinement 
energy is much greater than the interaction energy, then 
excitations to higher-order transverse modes are largely 
suppressed. All that the interaction can do in this case is 
multiply the lowest-order transverse mode by an almost 
uniform phase factor (a non-uniform phase would im- 
ply non-negligible amplitudes in higher-order transverse 
modes), thus allowing high-fidelity quantum gates. Suffi- 
ciently strong confinement could be achieved for example 
using hollow core photonic crystal fibers [l^ or optical 
nanofibers [soj . 



imposing important constraints on the performance of 
potential quantum gates. For contact interactions the 
effective phase is essentially always zero, no matter how 
strong the interaction. For long-range interactions the 
situation is more favorable, but there are still significant 
trade-offs between the achievable phase and fidelity. We 
discussed two potential solutions. One is to have a signif- 
icant transverse separation between the two wave pack- 
ets, which is possible only for long-range interactions. 
The price to pay is the need for an even stronger in- 
teraction. The second potential solution is to impose 
very strong transverse confinement, which may be possi- 
ble using hollow fibers or nanofibers. In any case it will 
be essential for future implementations of photon-photon 
gates to take transverse multi-mode effects into account. 



VI. SUMMARY 

The importance of the multi-mode character of quan- 
tum fields for the implementation of photon-photon gates 
had been recognized in the past, but only for the lon- 
gitudinal degree of freedom. Here we showed that the 
transverse degrees of freedom also play a significant role. 
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